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Abstract 

We discuss a scenario of the type-II 2HDM in which the 6677 rate of the Higgs pair production 
is enhanced due to the two nearly degenerate 125 GeV Higgs bosons {h, H). Considering various 
theoretical and experimental constraints, we figure out the allowed ranges of the trilinear couplings 
of these two Higgs bosons and calculate the signal rate of 6677 from the productions of Higgs pairs 
{hh, hH, HH) at the LHC. We find that in the allowed parameter space some trilinear Higgs 
couplings can be larger than the SM value by an order and the production rate of 6677 can be 
greatly enhanced. We also consider a ’’decoupling” benchmark point where the light CP-even Higgs 
has a SM-like cubic self-coupling while other trilinear couplings are very small. With a detailed 
simulation on the 6677 signal and backgrounds, we find that in such a ’’decoupling” scenario the 
hh and hH channels can jointly enhance the statistical significance to 5a at 14 TeV LHC with an 
integrated luminosity of 3000 fb“^. 

PACS numbers: 12.60.Fr, 14.80.Ec, 14.80.Bn 
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I. INTRODUCTION 


So far the properties of the 125 GeV Higgs boson discovered by the ATLAS and CMS 
collaborations jl, 2| agree with the Standard Model (SM) predictions. However, there is 
no experimental information for the Higgs self-conpling, which is vital for the spontaneous 
electroweak symmetry breaking. As is well known, the Higgs pair production at the LHC 
may provide a way to probe the Higgs self-coupling. The signal bbbb from the Higgs pair has 
the largest rate, but suffers from the huge QCD background. The bbrf channel is swamped 
by the bbjj background ^ where each light-flavored jet can fake a hadronic r. The detection 
of these two channels and also the bbWW* channel needs more elaborated strategies like 

Although the bh'j'j channel has a 




boosted kinematics and jet substructure technique 

small rate, it has the cleanest background, and thus has attracted more attention |5 Hio|. 
For the SM, the significance for gg ^ hh ^ 6677 is only around 2 ct at the 14 TeV LHC 
with an integrated luminosity of 3000 fb“^ js, 6 |. So a collider with higher energy (say 100 
TeV) seems needed to examine the Higgs self-coupling from the Higgs pair production. 

The Higgs pair production can serve as a good probe for new physics. The production 
rate can be enhanced by modifying the Higgs self-coupling or top quark Yukawa coupling 
properly. Also, it can be enhanced by some new mechanisms in the production, such as 
the heavy top partner loops in the little Higgs model [n|, the squark loops in the SUSY 
models [l2|, and the on-shell production of a heavy Higgs which decays into a pair of 125 
GeV Higgses in the two-Higgs-doublet model (2HDM) jsl, l9| . In this work, we will discuss a 
scenario in the type-H 2HDM |l^ where the bb'j'y channel of the Higgs pair is enhanced due 
to the two nearly degenerate 125 GeV Higgses (similar degenerate cases have been discussed 
in the literature [l^, but their impact on Higgs pair signals has not been studied). The mass 
splitting between these two Higgses is smaller than the mass resolution of the detector while 
larger than their widths so that the interference terms can be neglected. First, considering 
the theoretical constraints from vacuum stability, unitarity and perturbativity as well as the 
experimental constraints from the electroweak precision data, flavor observables and Higgs 
data, we will hgure out the allowed ranges of the trilinear couplings of these two Higgses 
and calculate the bb'j'j production rate at the LHC. Then, focusing on a ’’decoupling” 
benchmark point where the light CP-even Higgs has a SM-like cubic self-coupling while the 
other trilinear couplings are very small, we perform a detailed simulation on the bb'y'j signal 
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and its backgrounds at the 14 TeV LHC with an integrated luminosity of 3000 fb“^ . 

Our work is organized as follows. In Sec. II we recapitulate the type-II 2HDM. In Sec. 
Ill we describe our numerical calculations. In Sec. IV, we show the allowed ranges of the 
various trilinear couplings of the two Higgses and give the simulation results for the 6677 
signal and its backgrounds at the LHC. Finally, we draw our conclusion in Sec.V. 


II. TYPE-II 2HDM 


The general Higgs potential of 2HDM is written as 151 


V = + 17112 (^ 1 ^ 2 ) - mi 2 ($i*h 2 + h.c.) 




^($ 14 ) 2 ) 2 + h.c. 


+ 


A6($l<l>l)($I<l>2)+h.C. 


tr 


A 7 (<F|‘F 2 ) (‘^ 1 *^ 2 ) + h.c 


t, 


( 1 ) 


In the type-II 2HDM, a discrete Z 2 symmetry is introduced to make Ag = Ay = 0 while allow 
for a soft-breaking term with m ^2 7 ^ 0 - m ^2 taken to be real in order to avoid 

the explicit CP violation in the Higgs sector. 

The two complex scalar doublets have the hypercharge V = 1, 

V *.4 

(ui + 0? + lai) j y^(v2 + (j)l + ia2) 

where vi and V 2 are the vacuum expectation values (VEVs) with + n| = (246 GeV)^ 

and tau/d is dehned as v^jvx. The physical scalar spectrum of this model consists of two 
neutral CP-even h and iL, one neutral pseudoscalar H, and two charged scalar H^. This 
basis can be rotated to the Higgs basis by a mixing angle /9, where the VEV of $2 held is 
zero. In the Higgs basis, the mass eigenstates are obtained from 




h = sin(/9 — a)(f)^ + cos((3 — a) 02 ) 

H = cos(/5 — a)05 ~ sin(/9 — a)02, 

H = a 2 , = <Pt, (3) 


where the helds in the righ sides denote the interaction eigenstates in the Higgs basis. 
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as 


'n the Higgs basis, the general Yukawa interactions with no tree-level FCNC are written 
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Cy = — 



V 

where <hi(a;) 
have 


+ Kd^2)dR + MIQl{^i + Ku^2)ur + + K£^2)£r 


+ h.c., (4) 


= iT 2 ^*{x) and are the Yukawa matrices. For the type-II 2HDM, we 


Ku = cot 13, Kd = Ke = — tan (3. 


( 5 ) 


From Eq. (jl]) we can obtain the couplings of neutral Higgs bosons normalized to the SM 
Higgs boson 


Uy = sin(/3 — a), yj = sin(/5 — a) + cos{(3 — a)Kf, 
yy = cos(/9 — a), yf = cos(/5 — a) — sin(/3 — 

Vv = Vu = y^e = ( 6 ) 


where V denotes Z and W, and / denotes u, d and i. The charged Higgs couplings are 
written as 

Cy = VcKkiMdPR — Ku MuVckmPl] d + kMiPrI^ + h.c., (7) 

where Mj are the diagonal fermion mass matrices. 


III. NUMERICAL CALCULATIONS 


We employ 2HDMC-1.6.5 


17l | to implement the theoretical constraints from the vacuum 


stability, unitarity and coupling-constant perturbativity, and calculate the oblique parame¬ 
ters (S', T, U) and 5p. We use Superlso-3.4 [l^ to implement the constraints from B —)■ 


and use HiggsBounds-4.1.3 


19| to implement the exclusion constraints from the neutral and 


charged Higgses searches at the LEP, Tevatron and LHC at 95% conhdence level. The 
in-house code is used to calculate fit fo 125.5 GeV Higgs signal, Auib, and Aijib^- In 
addition to the theoretical constraints, we require the type-H 2HDM to satisfy all the exper¬ 
imental data at 2a level. The experimental values of electroweak precision data, B —)■ X^y, 
ArriBs Atub^ are taken from 201. 

We g enerate the 2HDM@NLO model using the tree-level 2HDM model and NLOCT pack¬ 


age 


2l|. The model contains the QCD R2 vertice and UV counterterms for the 2HDM, which 
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is based on the FeynRules 22| and UFO 23j frameworks. In our simulatiom 


signal and background events are generated with MadGraph5_aMC_v2.3.0 


he parton level 


24| . For the Higgs 


pair production via gluon-gluon fusion, we take the factorization and renormalization scales 
as the the invariant mass of the Higgs pair. The in-house code is used to transform the 
results of 2HDMC into the parameter card which is read by MadGraph5_aMC_v2.3.0 since 
there are different basis and mixing angles in CP-even Higgs sector between 2HDM@NLO 


model and 2HDMC. PYTHIA 


table of parameter card, and perform parton shower and 
fast detector simulations and data analysis with Delphes 


25l | is employed to decay the Higgs bosons following the decay 

ladronization. We perform the 


26 1 and Madanalysis5 27|. Jet 


reconstruction is done using the anti-fc-r algorithm with a radius parameter of i? = 0.5. The 
efficiency for 6-tagging is taken as 70%. The efficiency of photon tagging and the mis-tagging 
of QCD jets is assumed to the default value as in Delphes. 


Using the method in 




we perform a global £t to the 125.5 GeV Higgs data of 29 


channels after ICHEP 2014 29|]. Since we assume that the mass splitting of the two CP- 
even Higgses is smaller than the mass resolution of detector, the signal strength for a channel 
is dehned as 


E 


'-ggH^^ggH 

H=h, H 




( 8 ) 


where = (crx /(ax Hi?)with j denoting the partonic process FT, VBFH, VH, 
or ttH, and e* denotes the assumed signal composition of the partonic process j 30| , which 
has the same value for h and H. For an uncorrelated observable i, 


2 ^ (hi 
X-i 


A 


exp \2 


a/ 


( 9 ) 


where and a* denote the experimental central value and uncertainty for the i-channel. 
The uncertainty asymmetry is retained in our calculations. For the two correlated observ¬ 
ables, we take 


Y- = 


l-p 2 _ 


(ft-ft"")" , (ft 

o ~r 


A 


exp\2 

j 






-2p 


(A. - aD (a. - aD 


(^i 


a,- 


( 10 ) 


where p is the correlation coefficient. We sum over in the 29 channels, and pay particular 
attention to the surviving samples with — Xmin — 6.18, where Ymm denotes the minimum 
of x^- These samples correspond to the 95.4% conhdence level region in any two-dimension 
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plane of the model parameters when explaining the Higgs data (corresponding to the 2a 
range). 

In onr calcnlations, we take rrih = 125.5 GeV and rriH = 126 GeV, and the inpnt pa¬ 
rameters are cos(/5 — a), tan/3, the physical Higgs masses (m^, mH±) and the soft breaking 
parameter m^ 2 - Since the Higgs conplings between the two GP-even Higgses are independent 
of ruA and mH±, we take niA = which is favored by the 6 p and obliqne parameters. 

We scan randomly the parameters in the following ranges 

0 < cos(/3 — a) < 0.1 ,1 < tan /3 < 15, 

200 GeV <mA = mH± < 700 GeV, - (400 GeV)^ < ml^ < (400 GeVf. (11) 

We take the convention 0 < cos(/3 —a) < 1 and —1 < sin(/3 —a) < 1. With 0 < cos(/3 —a) < 
0.1, the conplings between the light GP-even Higgs and the gange bosons are close to the 
SM predictions while the corresponding heavy Higgs conplings are very small. 

IV. RESULTS AND DISCUSSIONS 

After imposing the above mentioned theoretical and experimental constraints, we hnd 
the minimal valne of is xlmin — 18-08, which is slightly larger than the SM valne (17.0). 
And the corresponding parameters are 

sin(/3 — a) ~ 0.99996, tan/3 ~ 3.094, mh = 125.5 GeV, mn — 126.0 GeV, 
mA = 448.88 GeV, mH± = 448.88 GeV, ml^ = 4615.4 GeV^. ( 12 ) 

A. Higgs pair cross section and Higgs trilinear couplings 

We dehne as the bb'yy signal event nnmber of type-H 2HDM normalized to the SM 
prediction 

Zcrjgg ^ HH) x BrjHH ^ 6677 ) .... 

a{gg hh)sM x Br{hh -)■ 6677 ) 51 ^’ 

where HH denotes hh, hH or HH. In fact, the contribntions from gg —)■ HH bb'y'j can 
be neglected since Br{H —)■ 77 ) is mnch smaller than the SM prediction for 0 < cos(/3 —a) < 

0 . 1 . 

In Fig. [H we project the snrviving samples on the planes of cos(/3 — a) versns tan /3 
and m ^2 versns tan/3, respectively. At the 14 TeV LHG with an integrated Inminosity of 
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FIG. 1: The scatter plots of surviving samples projected on the planes of cos(/3 — a) versus tan/3 
and m ^2 versus tan/3, respectively. The crosses (red) are for ^ 2^677 — 1-2) and bullets (green) for 
1.2 < <2.0, and triangles (black) for >2.0. 


3000 fb the significance of SM is aronnd 2a for the 6677 channel j^, . Therefore, it 

should be difficult to probe the 6677 channel of type-II 2HDM for < 2.0. As shown 

in this hgure, > 2.0 favors tan/3 < 2 (tan/3 < 1.2 is excluded by and Am^^), 

— lx 10^ GeV^ < m ^2 < —3 x 10^ GeV^ and 0 GeV^ < ml 2 < 1.5 x 10^ GeV^. The various 
Higgs trilinear couplings are sensitive to tan/3 and m^ 2 - addition, the top quark Yukawa 
couplings is sensitive to tan/3, as shown in Eq. ([ 6 ]). 

To understand the allowed ranges of we project the surviving samples on the 

planes of the Higgs couplings in Fig. [2j The upper panel of Fig. |2] shows that the light 
GP-even Higgs trilinear coupling and its coupling to top quark are restricted to be around 
the SM predictions, respectively. The absolute value of the heavy GP-even Higgs coupling 
to top quark is always suppressed, and allowed to be as low as 0.12 relative to the SM top 
quark Yukawa coupling. In some parameter space, the absolute value of the Higgs trilinear 
couplings of HHH and HHh are respectively allowed to be as high as 15 and 10 relative 
to the SM hhh coupling. The absolute value of the coupling Hhh is always suppressed 
compared to the SM hhh coupling due to the suppression of cos(/3 — a). 

From Fig. [2] we see that Rbi^^ > 2.0 favors two different regions. In one region, the Higgs 
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yHHh yHHh 


FIG. 2: Same as Fig. [H but projected on the planes of yf versus Hhhh, versus yHHH, Uhhh 
versus yHHh and ynhh versus ynnh- All these Higgs trilinear couplings are normalized to the SM 
hhh coupling. 

potential is ’’decoupling”, namely the hhh coupling is near the SM prediction while other 
trilinear couplings of HHH, HHh and Hhh are very small. Therefore, for the gg —)■ hH 
production process, the contributions of triangle diagrams will be sizably suppressed since 
the couplings of HHh and Hhh are very small. This will sizably soften the destructive 







FIG. 3: Same as Fig. [H but projected on the planes of Rhu versus Rhh and Rrh versus Rhh- Here 
the ratios denote the production rates via the gluon fusion normalized to the SM cross section of 
hh production. 

interference between the triangle and box diagrams, leading the enhancement of the cross 
section of gg hH. In the other region, the conpling HHh is mnch larger than the SM hhh 
coupling, which can make the contributions of the triangle diagrams to overcome the box 
diagrams, and enhance the cross section of gg —)■ hH. In addition, the upper-right panel of 
Fig. [2] shows that > 2.0 favors the absolute value of to be larger than 0.5, which 
avoids the cross section of gg —)■ hH to be sizably suppressed. 

Although the couplings oi HHH and HHh can be much larger than the SM hhh coupling, 
the cross section oi gg ^ HH can not be enhanced since there are destructive interference 
between the triangle diagrams mediated by H and h. Conversely, the cross section of 
gg —)■ HH is smaller than the SM cross section of gg —?• hh since the Htt coupling is 
suppressed. We show the cross sections of hh, hH and in Fig. [3l This hgure shows that 
the cross section of HH is smaller than 0.6 relative the SM hh prediction for most surviving 
samples. The cross section of hh is around the SM prediction, and the cross section of hH 
can reach 17 times of the SM hh prediction. 
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B. Simulation results in a decoupling scenario 


As seen from the preceding section, the cross section of Higgs pair production at the 
LHC can be sizably enhanced by a large Higgs trilinear coupling in the 2HDM with two 
nearly degenerate CP-even 125 GeV Higgs bosons, and as a result the Higgs pair signal is 
observable at the LHC. In the following we consider a ’’decoupling” scenario in which the 
light CP-even Higgs has a SM-like cubic self-coupling while other Higgs trilinear couplings of 
the two CP-even Higgses are very small. In this scenario, the bb —?■ hh, hH, HH processes 
can be neglected since there is no enhancement of Higgs trilinear couplings. 

We take a benchmark point 


sin(/3 - a) ~ -0.999988, tan/3 ~ 1.232, nih = 125.5 GeV, ~ 126.0 GeV, 
niA = 595.65 GeV, mH± = 595.65 GeV,m ^2 = 12304.0 GeV^; 
yhti = -0.996, ymt = 0.82, = -1.006, y^f^i = -1.228, 

Uhhh = —0.99996, yHHH = 0.245, ynnh = 0.0598, ynhh = 0.00552; 

Br{h 77) = 1.969 x 10■^ Br{h bb) = 0.6119, 

Br{H 77) = 9.702 x 10■^ Br{H bb) = 0.8385; 

a{gg -)■ hh)LO = 16.74 fb, a{gg hH)LO = 50.4 fb. (14) 

Since Br{H —)■ 77) is very small, we neglect gg —?■ HH —)• 6677, and consider the 6677 
signal from 


gg ^ hh ^ 6677 , 

gg —)■ hH —)■ 6677 . (15) 


The main SM backgrounds include non-resonant bb'j'j, tih {tt bb + X, h ^ 77 ), Zh 
{Z —)■ 55, h —)■ 77 ) and bbh {h 77 ). We neglect the subdominant reducible backgrounds 
of jj 77 and ff ^77 [ 8 |. The QGD corrections are considered by inclu ding a fc-factor, which is 
2.27 for the signal^M], 2.0 for bb'j'j jo], 1.1 for tth 32|, 1.33 for Zh 32| and 1.2 for bbh 33 1. 

Fig. m shows the distributions of some kinematical variables at the LHG with y/s = 14 
TeV for the hh, hH, bb'j'j and tth. The results oi Zh and bbh are not shown since they are 
subdominant. According to the distribution differences between the signal and backgrounds, 
we can improve the ratio of signal to backgrounds by making some kinematical cuts. First, 
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FIG. 4: The Higgs pair signal 6677 and background distributions of ^ Ai?b 6,77 and 
at the 14 TeV LHC. 


we require the final states to include two isolated photons and two 6-jets, and further impose 
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the following cuts 


> 60 GeV, > 25 GeV, P?' > 60 GeV, P]? > 25 GeV, 

> 0.4, > 0.4, > 0.4, 

\r]b\ < 2.5, \'q^\ < 2.5, 

Mbb > 30 GeV, > 30 GeV, Mbb^^ > 350 GeV, (16) 

where P!^^ and P^^ denote the transverse momentum of the hardest cind photon, and 
P!l? and P^^ for the second hardest &-jet and photon. Ai? = (A0)2 + (At])^ is the par¬ 

ticle separation with A(/) and A?] being the separation in the azimuthal angle and rapidity 
respectively. The cuts of the invariant mass of two 5-jets and two photon Mbirf-y, P^ and 
P]} can suppress the backgrounds sizably, especially the largest background 
The photon pair is further restricted to have 

/\R^^ < 2.0, 115 GeV < < 135 GeV. (17) 

The 5-quark pair is restricted to have 

/\Rbb < 2.0, 100 GeV < Mbb < 140 GeV. (18) 

Since the two photons (two 5 quarks) in the signals are from the Higgs decay, the signal 
rates peak at their invariant mass around the Higgs mass with relative small separation. 
The cuts in Eqs. ffT7)l and ffTSjl play the dominant role in suppressing the backgrounds. 

Finally, we make some cuts which can specially suppress the background tth. Since 
will decay into the charged leptons and jets, the background tih tends to include additional 
charged leptons and more jets. Therefore, we will veto the following case 

P^ > 20 GeV or > 20 GeV, (19) 

where Pf® denotes the transverse momentum of the 8th hardest jet. 

The resulting cut flow is shown in Table. [11 The bby'y and tth are the two major back¬ 
grounds. After imposing the above cuts, the events from signal hH are approximately 1.6 
times of those of hh. Since h and H are the degenerate 125.5 GeV Higgses, the total signal 
events are from hh and hH, whose signihcance can reach 5a at the 14 TeV LHG with an 
integrated luminosity of 3000 fb~^. If there is sizable mass splitting between h and H, the 
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TABLE I: The cut flow for the signal and background event numbers at the 14 TeV LHC with an 
integrated luminosity of 3000 fb“^ for the 6677 channel. The two culumns labeled ’/i/i’ and 'hH' 
are for the Higgs pair signal while other culumns are for the backgrounds. 


y/s = 14 TeV, 3 ab~^ 

hh 

hH 

6677 

tth 

Zh 

bbh 

s/Vb 

after cut in Eq.ffTHI) 

38.4 

63.3 

15999 

246.8 

22.5 

9.8 

0.8 

after cut in Eg.ffTTl) 

29.9 

48.1 

679.5 

152.5 

18.4 

4.7 

2.7 

after cut in Eq.ffTSl) 

18.7 

29.9 

74 

16.8 

2.7 

0.4 

5.1 

after cut in Ea.ffTOl) 

18.6 

29.7 

74 

10.4 

2.7 

0.4 

5.2 


cuts in Eq. ffTTl) and Eq. ffTSli will hurt the events of hH inevitably and hence suppress the 
signihcance. Therefore, the degeneracy between h and H plays the key role in enhancing 
the signihcance for such a ’’decoupling” scenario. 


V. CONCLUSION 

In this work we discussed a special scenario in the type-II 2HDM where the 6677 channel 
of the Higgs pair production can be enhanced due to the two nearly degenerate 125 GeV 
Higgses. We considered various theoretical and experimental constraints and found that 
in the allowed parameter space some trilinear Higgs couplings can be larger than the SM 
value by an order and the signal bh'y'y can be sizably enhanced. We also considered a 
’’decoupling” scenario where the light CP-even Higgs has the SM-like cubic self-coupling 
while other trilinear Higgs couplings are very small. From a detailed simulation on the 
signal 6677 and backgrounds, we found that the hh and hH production channels can jointly 
enhance the statistical signihcance to 5a at the 14 TeV LHC with an integrated luminosity of 
3000 fb“^. Therefore, the degenerate h and H play the vital role in enhancing the signihcance 
for probing the Higgs potential ’’decoupling” scenario. 
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